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Abstract: The expressions of solutions for general n x m matrix-valued inhomo- 
geneous linear stochastic differential equations are derived. This generalizes a result of 
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Jaschke (2003) for scalar inhomogeneous linear stochastic differential equations. As an 
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application, some M n vector-valued inhomogeneous nonlinear stochastic differential equa- 
tions are reduced to random differential equations, facilitating pathwise study of the so- 



1 A Review of Stochastic Exponential Formulas 

We first review some existing results about solution formulas for linear stochastic differen- 
tial equations (SDEs) or for their integral formulations. Let (£l,J-, {Ft), P) be a standard 
stochastic basis. For the following stochastic integral equation 

X t = l+ f Xs-dZ s , (1.1), 
Jo 
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where Z is a semimartingale with Z$ = 0, Doleans-Dade (1970) proved that the unique 
solution of (1.1) is given by 

X t = e W {z t -hz c ,Z c )t} [] (l + AZ s )e- AZ °. (1.2) 

0<s<t 

In the literature X is called Doleans exponential (or stochastic exponential) of Z, and 
is denoted by £(Z). The formula (1.2) is called the Doleans (or stochastic) exponential 
formula. 

In an unpublished paper, Yoeurp and Yor (1977) proved the following result for the 
solution formula of scalar SDEs (see also Revuz-Yor (1999) and Protter (2005) for the case 
where Z is a continuous semimartingale, and Melnikov-Shiryaev (1996)) for the general 
case). 

Theorem 1.1 (Yoeurp and Yor 1911 ) Let Z and H be semimartingales, and AZ S ^ 
— 1 for all s G [0, oo]. Then the unique solution of the inhomogeneous scalar linear SDE 

rt 



is given by 



where 



X t = H t + f X s _dZ s , t > 0, (1.3) 
J o 

X t = S(Z) t {H + J*S(Z)jldG s }, (1.4) 



G t = H t -{H^) t - £ ^^L. (1.5) 



Jaschke (2003) extended equation (1.3) to the case where H is an adapted cadlag 
process, not necessarily a semimartingale. He proved that in this case the solution of (1.3) 
is given by: 

X t = H t - £(Z) t f Hs-diEiZ)- 1 ). (1.6) 
Jo 

By using (1.6) Jaschke (2003) obtained a new proof of (1.4). 

On the other hand, Emery (1978) considered the following nxn matrix- valued stochas- 
tic equation 

U(t) = I+ f\dL(s))U(s-), (1.7), 
Jo 

where /isannxn identity, L is a given nxn matrix-valued semimartingale with L(0) = 0, 
such that / + AL(s) is invertible for all s G [0, oo], where AL(s) = L(s) — L(s-). Emery 
proved that the equation (1.7) admits a unique solution U, which isnxn matrix-valued 
semimartingale. We call it the stochastic exponential of L and denote it by £ (L). However, 
there is no explicit expression for such a stochastic exponential in general. 
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Jacod (1982) has studied the following inhomogeneous matrix-valued stochastic inte- 
gral equation 

X(t) = H(t)+ f (dL(s))X(s-), (1.8), 
Jo 

where L is a given n x n matrix-valued semimartingale with LO = 0, such that I + AL(s) 
is invertible for all s £ [0, oo], and H is an n x m matrix valued semimartingale. For an 
n x n matrix valued semimartingale A and an n x m matrix valued semimartingale B, we 
let 

[A,B](t) = (A c ,B c )(t) + AA(s)AB(s). 

0<s<t 

Here A c denotes its continuous martingale part defined componentwise by (L )^ = (L l j) c , 
and 

(A c , B c )) = ^2((A c yk, (B c )j). (1.9) 
k 

Using these notations the result of Jacod (1982) implies the following 
Theorem 1.2 (Jacod 1982) The unique solution of (1.8) is given by 

X(t) = S(L)(t){H(Q) + jy^is-r^Gis)}, (1.10) 

where 

G(t) = H(t) - (L c ,H c )(t) - (! + AL(s))- 1 AL(s)A J ff(s). (1.11) 

0<s<t 

In particular, if L and i? are continuous semimartingales, an expression for the solution 
of (1.8) is given in Revuz-Yor (1999) as follows: 

X(t) = S(L)(t){H(0) + J* SiL^sr^dHis) - d[L,H](s))}. (1.12) 

The objective of the present note is to generalize equation (1.8) to the case where H{t) 
is a given n x m matrix-valued adapted cadlag process, not necessarily a semimartingale. 
We give an expression of the solution of (1.8) for this case and also give a simpler proof 
for Theorem 1.1. Our result extends (1.6) of Jaschke (2003) to matrix-valued case. As 
an application, we reduce some iR ra -valued inhomogeneous nonlinear SDEs to random 
differential equations (RDEs) — differential equations with random coefficients. This 
facilitates pathwise study of solutions and is an important step in the context of random 
dynamical systems approaches; see Arnold (1998). 
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2 Matrix-valued Inhomogeneous Linear SDEs 

Leandre (1985) obtained the following result about stochastic equation (1.7). If we denote 
by V the inverse of S(L), then V is the solution of the following equation: 

V(t) = I+ f V{s-)dW{s), 
Jo 

where 

W(t) = -L(t) + (L c ,L c )(t)+ ]T (1 + AL( S ))^(AL( S )) 2 . 

0<s<t 

That means £{L)~ l = £(W T ) T . We refer the reader to Karandikar (1991) for a detailed 
proof of this result. 

Now we will use this result of Leandre (1985) to solve the following inhomogeneous 
stochastic integral equation 

X(t)=H(t)+ f\dL(s))X(s-), (2.1), 
Jo 

where L is a given n x n matrix- valued semimartingale with LO = 0, such that I + AL(s) 
is invertible for all s € [0, oo], H(t) is a given nxm matrix- valued adapted cadlag process. 

Our main result is the following. 
Theorem 2.1 The unique solution of (2.1) is given by 

X{t) = H{t) - £{L){t) f (d£(L){ S y 1 )H{s-). (2.2) 

Jo 

If H is n x m matrix-valued semimartingale, then X(t) has the same expression as given 
by (1.10) and (1.11). 

Proof. We denote £(L) and £(L) -1 by U and V, respectively. We are going to show 
that the process 

X(t) = H(t) - U{t) f\dV( S ))H(s-), 
Jo 

defined by (2.2), satisfies equation (2.1). Since UV = I, by the integration by parts 
formula (see Karandikar (1991)) we get 

= d(U{t)V{t)) = U{t-)dV(t) + (dU{t))V{t-) + d[U, V](t) 
= d{U{t)V{t)) = U{t-)dV{t) + dL{t) + d[U, V](t). 

Once again by the integration by parts formula, using the above result and the fact that 
dU{t) = (dL(t))U(t-), we have 

d(X(t)-H(t)) = -U(t-)dV(t)H(t-) - (dU(t))( C dV(s)H(s-)) - (d[U,V](t))H(t-) 

Jo 

= (dL(t))[H(t-) - U(t-)( dV(s)H(s-))} = (dL(t))X(t-). 

Jo 
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This shows that the process (X(t)) defined by (2.2) satisfies (2.1). 

Now we assume that (H(t)) is an n x m matrix-valued semimartingale. Using the 
notations in Section 1 we can verify that 

G(t) = H(t) + [W,H](t). (2.3) 

By the integration by parts formula, using (2.3) and the fact that dV(t) = V(t—)dW(t), 
we obtain 

= d(V(t)H(t)) = V(t-)dH(t) + (dV(t))H(t-) + V(t-)d[W, H](t) 
= V(t-)dG(t) + dL{t) + (dV(t))H{t-), 

from which we get 

H(t) = U(t){H(0) + £v(s-)dG(s) + J* dV(s)H(s-)}. 

Thus, if we let (G(t)) be defined by (2.3) then (X(t)) has the expression of (2.2), and 
consequently it satisfies equation (2.1). The proof of the theorem is complete. □ 



3 An Application to Nonlinear SDEs 

In this section we will apply our results in Theorem 2.1 to reduce an inhomogeneous 
nonlinear SDE to a RDE (random differential equation). Now we consider the following 
n-dimensional nonlinear SDE (but with a linear multiplicative noise term): 

n 

dX\t) = f(t,X(t))dt + Y,Cj(t)X j (t)dB j (t), X i (0)=x i , (3.1) 

where C(t) is an n x n matrix-valued measurable function, f(t,x) is a iR ra -valued measur- 
able function on [0, oo) x M n , and B(t) = {B 1 ^), B n (t)) T is a n-dimensional Brownian 
motion. Put 

L){t) = J* C^dB^s), i.j = l,---,n; H(t) = X(0) + jf* f(s,X(s))ds. (3.2) 

Then (3.1) can be rewritten in the form of (2.1). For such L, the equation (1.9) is reduced 
to the following linear equation: 

n 

dU l j (t) = J2Cl(t)U*(t)dB k (t), Up) = 5% (3.3) 
k=i 

In the present case we have G(t) = H(t). So according to (2.2), the solution of (3.1) can 
be expressed as 

X(t) = U{t){x + jf* Uit)- 1 /^, X(s))ds}, (3.4) 
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where U is the solution of (3.3). Unfortunately, even in this case we are not able to give 
an explicit expression for U(t). Let Y(t) = L r (t) _1 X(t), then 



This is the integral formulation of a RDE (random differential equation). Once we have 
sample path solution Y(t) of this transformed RDE , we obtain the solution of the original 
SDE via X(t) = U(t)Y(t). 
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